Plasma dispersion of multisubband electron systems over liquid 

helium 

Sviatoslav S. Sokolov 

Departamento de Fisica, Universidade Federal de Sao Carlos, 13565-905 Sao Carlos, Sao Paulo, 
Brazil and B. I. Verkin Institute for Low Temperature Physics and Engineering, National 
Academy of Sciences of Ukraine, 61164 Kharkov, Ukraine 

Nelson Studart 

Departamento de Fisica, Universidade Federal de Sao Carlos, 13565-905 Sao Carlos, Sao Paulo, 

Brazil 

Density-density response functions are evaluated for nondegenerate multisub- 
band electron systems in the random-phase approximation for arbitrary wave 
number and subband index. We consider both quasi-two-dimensional and quasi- 
one-dimensional systems for electrons confined to the surface of liquid helium. 
The dispersion relations of longitudinal intrasubband and transverse intersub- 
band modes are calculated at low temperatures and for long wavelengths. We 
discuss the effects of screening and two-subband occupancy on the plasmon spec- 
trum. The characteristic absorption edge of the intersubband modes is shifted 
relatively to the single-particle intersubband separation and the depolarization 
shift correction can be significant at high electron densities. 

PACS number(s): 73.20.Dx; 73.20.Mf; 73.90.+f 

I. INTRODUCTION 

Two-dimensional electron systems (2DES) over the liquid helium surface have been in- 
tensively studied for a long time. More recently, it was possible to confine these surface 
electrons (SE) in reduced geometries creating also one-dimensional (ID) systems. Both 
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systems provide a nearly ideal laboratory for studying collective phenomena in the electron 
plasma in lower dimensions because the cleanness of the helium surface restricts SE scat- 
tering mechanisms to those with helium atoms in the vapor phase, which predominates at 
T > 1 K, and with surface oscillations (ripplons) at lower temperatures. Furthermore both 
scattering mechanisms become ineffective with lowering temperature and can be discarded at 
T < 0.1 K. In such a regime, collective effects in low- dimensional electron systems due to the 
Coulomb interaction can be investigated ignoring the interaction with scatterers. Another 
important feature of these systems is the accessible range of SE densities which is limited 
to Us < 10^ cm~^ (for bulk helium). As a consequence, the 2D Fermi energy Sp 10"^ K 
and SE behave like nondegenerate low-dimensional systems differing in many aspects from 
its quantum counterpart realized in semiconductor structures. 

As it is well known, the collective excitation spectrum depends crucially on the way the 
particles are confined. For instance, for longitudinal plasma oscillations of the 2DES, the 
spectrum uj2d{(1) ~ 0^^"^ is in contrast to the 3D situation in which one has a optical mode 
starting from the plasma frequency. This is a consequence of the fact that the screening is 
incomplete in 2D because there are electromagnetic fields in the vacuum surrounding the 
plane and many-body effects play important role in describing the properties of the 2DES. 
On the other hand, the longitudinal plasmon mode in the IDES case is oJiD^qx) ~ <ix ln(ga;^), 
where £ is some characteristic length of the system. In these cases, we have assumed that only 
the lowest subband, for electron motion along the direction perpendicular to the electron 
sheet, is occupied. This limiting case is achieved when the Bolztmann factor exp(— A21/T) ^ 
1, where A21 = A2 — Ai is the energy gap between the lowest (1) and the first-occupied (2) 
subband, and the occupation of higher subbands is negligible. Otherwise, the multisubband 
nature of low-dimensional electron systems — hereafter referred as quasi-2D(lD)ES — cannot 
be discarded when the temperature is comparable with A21 and population effects of higher 
subbands cannot be ignored. 

In this paper, we address the problem of plasmon spectrum in Q2DES and QIDES 
over the surface of liquid helium. We use the many-body dielectric formalism within the 
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random-phase approximation. In this approach, the mode spectrum is obtained from the 
roots of a determinantal equation for the dielectric function. At first glance, we note that 
the multisubband character of these systems allows the existence of transverse modes of 
plasma oscillations in the direction normal to that of unconfined electron motion. 

We adopt a two-subband model where the bare electron-electron interaction is evaluated 
using subband wave functions found by the variational method for the Q2DES and taken 
as the harmonic-oscillator functions for the parabolic confinement in the QIDES. We limit 
ourselves to the case of low enough temperatures which allows us to disregard the coupling 
of plasma oscillations with ripplon modes. We do not also consider the possible transition of 
the electron system to the ordered state where the electron-ripplon interaction can strongly 
modify the mode spectrum. [^|-|^ 



The main theoretical approach to the study of plasma oscillations in multisubband low- 
dimensional charge system is based on many-body dielectric formalism using the generalized 
dielectric function 



where Umm'i^, q) is the density-density response function, 5„„/ is the Kronecker symbol and 
Vnn'.mm'{(l) IS the matrix element of Fourier-transformed Coulomb interaction averaged over 
wave functions of subbands with indices n, n', m, and m' equal to 1, 2, 3.... The dielectric 
function enn^mm'iy^i q) depends both on the frequency and the wave numbers q for the Q2DES 
and qx for the QIDES. 

In the random-phase approximation (RPA), we assume that the electron system responds 
to external perturbations as a noninteracting system and we take limm'i!^, l) = nj^„, (u;, q), 
where the free polarizability function is written as 



II. THEORETICAL APPROACH 




(1) 
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Here i^k = h'^k'^/2m, where m is the electron mass, S is infinitesimal positive, and a is spin 
index. For classical systems, the distribution function /o(-Ek + ^n) = exp [— + A„) /T] 
and is normalized by the condition Xln k o- /o(-£'k + A„) = where N is the number of 
particles. 

Using the dielectric function given by Eq. (|T|), Vinter and Das Sarma |^ have studied 
many-body effects in the degenerate Q2DES. Das Sarma and co-workers ||9HTT|, Hu and 
O'Connell [l^ and Hai et al. [|l3l extended these studies to plasma oscillations in degenerate 
QID multisubband system whereas Sokolov and Studart fl^ approach the problem in the 
classical regime. 

The well-known bare electron-electron potential is given by 



f*oo poo 

VnS,mm'{Q) = I / dzdZ'1pniz)i)n' iz)v^" {q)lp^{z')ij^, (z) , (3) 
-'0 "'0 



and 



r'OO POO 

ynn'mm'iQx) = I / dydy' ^n{y)^n'{y)v^^ {qx)^m{y')Vm'{y') , (4) 



— oo J — oo 



where f^^(g) = 2Txl?/Sq, \v^^{qx) = 2{1? / L^)KQ{\q^\\y—y'\) is the Coulomb potential, S [L^] 
the area [length] of the system, and ipn{z) [^Pn{y)] denotes the n-th subband wave functions 
for the Q2DES [QIDES]. Here e = [2e^/(l + with e the helium dielectric constant, is 
the effective charge taking substrate effects into account. 



III. PLASMON SPECTRUM 

The dispersion relations for collective modes for a multisubband system are found from 
the roots of the determinantal equation 

det \enn', mm' {q, ^)\ = 0- (5) 

In principle, all the subbands should be considered in the above equation. However an useful 
analytical solution is possible in a two-subband model. In this case, Eq. (^) splits into two 
independent equations 
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(6a) 



1 - Vi2,i2 [n?2(^, q) + (a;, g)] = 0. (6b) 

Mode coupling appears only if one take into account higher subbands. Equation (^) de- 
scribes the longitudinal intrasuhhand plasma oscillations whose dispersion law must coincide 
with that of 2DES or IDES with one-subband occupancy system whereas Eq. ( |6b| ) gives the 
dispersion law for transverse intersubband oscillations involving transitions from the lowest 
to the second subband. 



A. Q2DES 



As it is well-known, SE on helium are trapped in the direction perpendicular to the 
surface {z direction) by a potential well due to image forces and a holding electric field E±. 
For = 0, and the image potential V{z) = —Aq/z, where Aq = (e^/4)(£: — l)/{e + 1), and 
infinite potential barrier at the interface, the solution of the Schrodinger equation is given 
by [1^13 



n 3/2 



n 



5/2 



exp 



fioz\ (1) f 2koz 



n 



n 



(7) 



where kq = mKo/h^ (= 3/(2 (^)o), where {z)^^ is the mean electron distance from the plane) 
and L!"n\x) are the associated Laguerre polynomials. The energy subband is given by the 
hydrogen-like spectrum A„ = Ao/n^ where Aq = h'^KQ/2m. If the pressing electric field 
E± is turned on, there is no general analytical solution and we assume trial wave functions 
corresponding to two lowest subbands (n = 1 and 2) of Eq. (|^ with variational parameters 
Ki and K,2'- ||T8| , p!9 



?/'i(z) = 2K^^^2;exp {—Kiz) 



(8a) 



5/2 



/«12 



(8b) 



and subband energies: 



Ai = ^ - Ao/«i + -7: : 



(8c) 



A, 



Here k? 



6m 



1 + 



K 
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1 + 



2k,1 — fi;ifi;2 



12 



+ 



2^ 



1 + 



4k,1 — + ^2 
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(8d) 
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K,iK,2 + /^i- If W6 define ki(-E'±) = "'71^0 and k,2{E±) = ri2Ho, one can find t]i 



and ri2 as the roots of the system of equations given by 



^1 =0, 

to 



(9a) 



^2(^1 - '^ViV2 + ^57]^V2 - 11^1^2 + 7%) - 2^2(^1 - ^ViV2 + 10?7ir/2 - + 



— ) {5r]t - 10r/3r/2 + ISr^^r/^ - Ar]ir]l + 2^) 







(9b) 



where k,± = {3meE±/2h?Y^^ . For E± = 0, Eqs. (0) and ( pb]) reproduce the results, given 
by Eq. (|^) and respective eigenenergies with rji = 1 and 772 = 0.5. The numerical values of 
the variational parameters as a function of the pressing electric field are plotted in Fig. 1. 
We observe a rapid increase at low field and an asymptotic linear behavior at larger fields. 
With these values for r/i and 772? we depicted in Fig. 2, the field dependence of the energy 
of the lowest-subband and the energy gap A21. 

Using Eqs. (^, ( paj ) and ( pb| ) one can calculate the values of Vu^u and Vi2,i2 up to 
second-order in the parameters q/ni <til and q/{Ki + K2) ^ 1 as 



1,11 



1 _ ^ . 

4ki 4fi:2j ' 



(10a) 



V^2,i2 = v'^{q)a{E^)^ 



16g 



+ 



Iq" 



(10b) 



5(ki + K2) (ki + /t2)2_ 

with a{Ei_) = 60r]fr]2/[{rii + ri2y{rH — 771772 + rj^)]- The well-behaved form of a{E±^) does not 
influence strongly Vi2,i2 because a{E^ = 0) = 0.146 and a{Ei_) increases by increasing E^ 
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until reaches a maximum Omax = 0.281 near E±_ = 0.3 kV/cm and slightly decreases to 0.227 
at E± = 3 kV/cm. 

For the Q2DES, the noninteracting density-density response function, Eq.(^, can be 
calculated in a straightforward way. The result is 

N 



hquTZn 



exp i-AjT) U { CU ] - exp (-A„VT) U [ CJ 



■(+) 

■>nn' 



(11) 



where (l^} = [u + (A„ — A,„/) /h] /qux ± hq/2'muT, with ut = \/lT jm being the thermal 
velocity and = Xln^^P {.~An/T). Similar general structure of 11°^, (cj, q) is found in the 



classical regime of the electron gas in 3D case. The function U {Q is given by the integral 
1 r"^ o^T^('_„2^ /-C 



f/(C) = ^ / ^^^PLZZ = _2 exp(-C^) / exp(t2)dt + zv^exp(-C') 



(12) 



For n = n' and for small g ^ (2muj/hY^'^ Eq. ([TT|) can be approximately expressed through 
W{(), the well-known function in the plasma theory, as pT|-|2^ 



n;°!.(i^,9) 



.JLw (^) exp (-^ 

TZn \quT J \ T 



(13) 



The function W{C,) is connected with f/(C) by the relation 

W{C) = -{dU/dO/2 = l + CU{C)- 

Putting oj = uj{q) — i'jq, and assuming uj{q)/quT ^ 1 and \uj{q) — /S.21 j h\ j quT ^ 1, 
we obtain, in the two-subband model, the dispersion relation for longitudinal intrasubband 
mode frequencies (a;;) and damping (7^): 



\ ^ 3g 3g 
k£, 4ki 



(14a) 



TT- -exp 

{qury 



{qurY 



(14b) 



where (^2d{q) — {2e^/ma?')q and ko = 2e^/Ta^ is the 2D Debye wave number and a = 
(7™^)^^/^ is the mean interelectron spacing. One can see, from Eq. ( |14a| ) , that the first 
two-terms of the real part of the longitudinal branch are the same as in the classical 2DES 



2T|^3[. The extra term (3g/4/ti) comes from the effect of the layer thickness, because ki is 
related to the distance of the electron in the lowest subband from the plane, and lowers the 
dispersion relation at small wavelengths. From Eq. ( |14a| ), we also observe the competition 
between screening and layer thickness effects on the Q2DES. The Landau damping 7z(g) of 
the plasma oscillation is practically the same as in the 2DES. 

For transverse intersubband modes, we obtain the frequencies (ut) and damping (74): 



(q) - ^21 + 2t^2lt^2D 



(sh) 



(1 + 34^^V2^2i) Tg2 
1 + ^ : 



m uj. 



(sh) 
■'2D 



Ki + K2 



it{q) 



where UJ21 = A2i/^ and cu, 



[uJt{q) - UJ21] 
TT — exp 



quT 



[uJt{q) - UJ21] 



(15a) 



(15b) 



(sh) 
2D 



2a{E / hKQa . As one can see the transverse plasmon 



mode spectrum, given by Eq. (|15aj) , is quite different of the longitudinal branch and has a 
gap at g = 0. The frequency of the characteristic absorption edge is shifted, relatively to the 
frequency ljJi2 of the intersubband transition, shown in Fig. 2, by Aa; = \J uj^^ + 2u2iIjJ2]j — 
a;2i, which is the manifestation of the depolarization shift effect in transverse oscillations of 



the many-body system. |12] The experimental observation of Atu, should be very interesting 

by evidencing the role of Coulomb effects in the collective electron motion along the z 
direction. Note that for the 2D plasma parameter F = e^/aT < 1, where RPA is formally 
valid, (^20 ^21 and Acj ~ 0^2^ i.e. the absorption edge is very close to UJ21 being only 
slightly shifted to higher frequencies. By increasing g, ujt{q) decreases according the last 
term in brackets in Eq. (|15aj) . However our estimates show that, for T ~ 0.1 — 1.0 K 
and g ~ 10 — 10^ cm^^, the coefficient of the quadratic term is larger than that of linear 
one. However, in the long wavelength limit, these coefficients are so small that ujt{q) ~ 



"^21 Y 1 + '^^2D l^'i'L- longitudinal mode, 7t(g), given by Eq. ( |15b| ), is exponentially 

small such that |7t(g)| <^ '-^t(g)- 

The absorption edge of the mode ujt{q) depends strongly on co'21. For very small E'x, 



the experimental values of UJ21 are close to 3Ao/4/l and increase linearly with Ej^. For 



arbitrary Ej_, uj2\ is obtained from the gap energy, displayed in Fig. 2, and is in the range 
of 100 GHz to 1 THz. The polarization shift Au; ~ ^fh^ ~ "^s for T < 1? even though 
lAco"! ^ For example for rig = 10^ cm~^ and E± = 0, we estimate oo^^ff^ ~ 100 MHz 
^ UJ21 ■ This makes very difficult the direct experimental observation of depolarization shift 
at this electron density. However, the effect should be observable at higher densities (for 
instance Ug ~ 10^ cm~^) even though our results can not be quite reliable in this regime since 
RPA should not be valid in such a strongly correlated classical Q2DES. However, one can 
hope that the nature of plasma oscillations does not change drastically, at least qualitatively, 
even in the high density limit (F > 1) and the depolarization shift should be measured for 
the Q2DES on the helium surface. 



B. QIDES 

We now consider plasma oscillations in the QIDES created along a channel filled with 
liquid helium. As in previous work, [0,^] we consider a parabolic confinement U{y) = 
mujQy'^/2 with the frequency coq = {eE±mRy^'^, where R is the curvature radius of the 
liquid in the channel. Typical values of R vary from 10~^ to 10~^ cm. The spectrum 
of electron subbands along the y axis is En = tvjjQin — 1/2), n = 1, 2, 3... in addition to the 
subbands along the z direction. The motion along the x direction (the channel axis) is free. 
The frequency ujq increases with E_[_ achieving 100 GHz at E^ = 3 kV/cm for R = 5 x 10~^ 
cm. As ujq <^ UJ21, the multisubband system in transverse directions can be decoupled and 
we ignore electron transitions in z direction which are the same as discussed above. 



The noninteracting density- density response function was calculated in Ref. |jTj]. The 
result was 

2N 



exp [-{n - l)nuJo/T] U (cU) - exp [-{W - l)nuJo/T] U {^C' 



(16) 



hqxUT [1 + coth (^o/2T)] 

where c!^i = {uj /q^Ux) [1 + {ujQ/uj){n — n')]±hqx/2muT- For Hloq ^ T, when only the lowest 
subband (n = 1) is occupied, the expression for the response function is greatly simphfied 
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yielding no„,(cu,g,.) = -(iV/T)W(u;/g,MT). 

Using the wave functions of the two lowest subbands {n = 1 and n = 2) 



My) 



1 



exp 



My) 



V2 I r 

-yexp 



1/4 1/2 \ 27/2 / ' 3/2^"^^^^ \ 27/2 / ' 

where = (/i/mcuo)^''^, we obtain the matrix elements of Coulomb interaction from Eq. 
(I): 

1 



and 



^11,11 fe) = exp ^ ^ 



In- 



^12^2(9.) 



1_ 

2Lrf 



exp 



<ilyl 



<ilyl 



\<ixyQ\ 



for |ga;?/o| < 1 



qxyo\ 



(17a) 



(17b) 



52 



for |ga;2/o| < 1- 



2 \qxyQ\ 
where Wa,i3{x) is Whittacker function. 

Using Eqs. (H), (|6aD , and (|17a|) , taking oj = uj{qx) — i'Jq, and assuming uj{qx) / qxUx ^ 1 
and \uj{qx)—ujQ\/qxUT ^ 1 we obtain the dispersion relation of the longitudinal intrasubband 
modes in the long wavelength limit, \qxyf)\ <^ 1, as 



'^Mx) = f^ln- 



mi IqxVo 



iMx) 



exp 



qlyl 
4 



■ 3T£, _i 1 



qxyo\ 



^t{qx) 

{qxUrf 



exp 



^t{qx 

(qur) 



(18a) 



(18b) 



where 



n, 



(N/Lx) ^ is the mean interelectron distance along the channel. 



The longitudinal spectrum mode, given by Eq. ( |18a| ), has the same structure of the 



obtained previously in Ref. ||T4| and in Refs. I^Q],^ where a quasi-crystalline approximation 
was employed. However we found an additional second term in brackets, which should 
be quite small for reasonable values of T and i. Note also that condition uj/qxUT ^ 1 
assumed here is equivalent to T ^ e^/£ in the quasicrystalline approximation. It worth 
emphasizes that the present result was obtained within RPA which is valid in the opposite 
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limit T ^ e'^/i. Our conclusion is that the plasmon spectrum in the classical QIDES has 
little dependence on the plasma parameter and RPA results should be probably correct in 
wide range of electron densities. 

The transverse branch of collective excitations is rather interesting. Following the same 
steps as before, we arrive to 

1 



2 ^^l 



m£ 



In 



1 9x1/0 1 



(19a) 



1 + 



(l + 3a;;5V2cc'o) T (l + ^[JV^o) h 



m uj 



(sh) 

ho 



2muj 



(sh) 
ID 



TT- 



QxUt 



exp 



(qxUt) 



(19b) 



Here uj['^^^ ='e^/M. The first two terms in Eq. ( |19a] ) correspond to the result obtained in the 
quasicrystalline approximation PU|PU| if i/q is replaced by i in the logarithmic factor. The 
next term is the depolarization shift correction increasing the absorption edge frequency by 
Aio = \J uJq + 2ujquj^^I^^ — ujq ~ oJ^il^^ when oJ^ij^^ <^ oJq- One can see that the instability of the 
transverse mode {ujl{qx) < 0) in the limit of zero confinement {uoq = 0) is still manifested in 
our treatment. We call the attention that we found a quite different result in our previous 
work [jl4[ because we used an approximate expression Vl,^^2ilx) —'^/Lx considered in Ref. 
|T^. One estimative is that the polarization shift correction should be quite small for 
rii ~ 10^ - 10^ cm-i and T ^ 10"^ - 1 K such that e'^/i < T. For instance, c^SJ^ ~ 10 GHz 
for = 3 kV/cm and ni = 10^ cm~^ whereas Uq = 100 GHz at i? = 5 x 10~^ cm. However, 
this density range can not be achieved in experimental conditions. For higher densities, Aa; 
should be of the same order of c^o and the polarization shift should be observed even our 
results are based on the RPA. 
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IV. CONCLUDING REMARKS 



In the present work, we have used the many-body dielectric formahsm to calculate the 
spectrum of plasma oscillations for the classical Q2DES and QIDES formed on the liquid 
helium surface. We have obtained the general expression for the density-density Q2D and 
QID response functions for any frequency and wave number within the RPA. The results 
are valid at low temperatures since we have used a two-subband model in which only the 
lowest subbands of the motion in the direction normal to the electron layer (Q2D) and of the 
motion in direction across the conducting channel (QID) are occupied. The plasma disper- 
sion relations were found from the zeros of the determinantal equation for the generalized 
multisubband dielectric functions. We have obtained corrections to the gapless longitudinal 
modes, beyond the q^/^-behavior in the Q2DES and sound-like behavior, within logarith- 
mic accuracy, in the QIDES. The intersubband transverse collective frequency is higher 
than the corresponding single-particle excitation frequency both in Q2DES and QIDES. 
The absorption edge frequencies are increased by the depolarization shift which can be large 
at high densities. In this connection, the experimental study of intersubband transition 
in low-dimensional electron systems over liquid helium seems to be attractive, because of 
the accessibility of wide range of charge concentrations and low temperatures, to observe 
collective effects on spectroscopic transitions. PBI 



We conclude by pointing out some limitations of our approach. The results are based on 
the RPA, which works quite well at small values of the plasma parameter. We know that 
RPA results become worse as the dimensionality is reduced, but we do not know how to go 
beyond RPA in a controlled way mainly in the QIDES. We are, however, encouraged by 
the good agreement of our RPA results for the mode spectrum and those obtained in the 
quasi-crystalline approximation that is valid in the opposite limit of high plasma parameter. 
Other fact is the excellent agreement between the RPA theory and experiment on collective 



excitations in semiconductor quantum wells and wires ||3^. Our use of a two-subband 
model can be and should be improved in more realistic calculations p^. But we do not 
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expect the correction of including other subbands to be quahtatively significant though at 
low temperatures. 
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FIGURE CAPTIONS 



Fig. 1. Variational parameters r/i (straight line) and 7]2 (dashed line) as a function of the 
pressing electric field E± evaluated numerically from Eqs. (pa]) and (pE|). 

Fig. 2. Lowest-subband energy (straight line) and energy gap A21 (dashed line) of single- 
electron spectroscopic excitation of the Q2DES as a function of the pressing electric field 
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Fig. 1. 

Sokolov and Studart, 
"Plasma dispersion..." 




Fig. 2. 

Sokolov and Studart, 
"Plasma dispersion..." 



